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Abstract. We perform a renormalisation group analysis of the SO(IO) model proposed by Chang, Masiero 
and Murayama, which links the large atmospheric neutrino mixing angle to loop-induced transitions be- 
tween right-handed 6 and s quarks. We compute the impact on Bg — Bs mixing and find that the mass 
difference in the Bs system can exceed its Standard Model value by a factor of 16. 

PACS. 12.60.Jv supersymmetric models - 12.15.Ff Quark and lepton masses and mixing 
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1 Introduction 

Grand Unified Theories combine quarks and leptons into 
symmetry multiplets. This opens the possibihty to find 
links between the flavour structures of both sectors. Chang, 
Masiero and Murayama (CMM) have proposed an inter- 
esting supersymmetric SO(IO) model in which the large 
i^fi~t^T mixing angle can affect transitions between right- 
handed b and s quarks through supersymmetric loop dia- 
grams The SO(10)-symmetric superpotential has the 
form 
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Here 16 is the usual spinor comprising the matter su- 
perfields and the other fields are Higgs superfields in the 
indicated representations. The Yukawa coupling is a 
symmetric 3x3 matrix in generation space containing the 
large top Yukawa coupling. The two dimcnsion-5 terms in- 
volve the Planck mass Mpi and further Higgs fields in the 
indicated representations. The last term in (Q) generates 
small neutrino masses via the standard see-saw mecha- 
nism, once 16h acquires an SO(10)-breaking VEV. From 
jjit S> rric we observe a large hierarchy in Y^ , which must 
be largely compensated in Y^ in order to explain the ob- 
served pattern of neutrino masses. This is achieved in a 
natural way by invoking flavour symmetries, the simplest 
of which render Y^ and Y^'^ simultaneously diagonal. In 
this basis the remaining Yukawa matrix Y^ has the form 
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Here Vckm and J/pmns are the CKM and PMNS matrices 
encoding flavour mixing in the quark and lepton sectors 
pi and certain diagonal phase matrices have been omitted 
in (O. From (0) one realizes that the flavour structure 
of Y^ is neither symmetric nor anti-symmetric. This is 
possible, because the corresponding dimension-5 term in 
transforms reducibly under SO(IO). At some scale Mio 
between the Planck and GUT scales the 45h and 16h 
acquire VEVs and SO(IO) is broken to SU(5). The SU(5) 
superpotential reads 
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16 
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PI 



(3) 



The matter supermultiplets <P and N are the usual 10, 
5 and 1 from the decomposition of the 16. Furthermore, 
Y^ oc y^W45/Afpi and W45 and ujg are the VEVs of the 
45h and 16h fields. The SO(IO) Higgs fields comprise the 
SU(5) ones as IOh 13 5h and 10^ D 5h- The remaining 
components of the SO (10) Higgs fields are assumed heavy 
with zero VEVs. Finally, at the GUT scale the SU(5) the- 
ory is broken to the MSSM and the MSSM Higgs fields 
Hu C 5h and Hd C 5h couple to up- and down-type 
fermions, respectively. 

The soft SUSY-breaking terms are assumed universal 
near the Planck scale. The large Yukawa coupling in Y^ 
now renormalises the squark mass matrix. The renormal- 
isation group (RG) flow down to Afio (and to AfcuT) 
will keep its diagonal form (in the basis in which Y^ 
and are diagonal), but will split the mass of the 
third from those of the first two generations. The diag- 
onalisation of Y^ involves the rotation of in Q with 
f^PMNS- Since <!> unifies the left-handed (s)leptons with 
righthanded (s)quarks, the large atmospheric neutrino mix- 
ing angle will appear in the mixing of right-handed b and 
s squarks. 

We have performed a complete RG analysis of the 
CMM model and computed its impact on Bg — Bg mix- 
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ing. The key parameter in the analysis is the top Yukawa 
couphng yt, which drives the bR-SR mixing effect. The 
fact that the RG evolution of this coupling is governed 
by infrared quasi-fixed points in SO(IO), SU(5) and the 
MSSM allows us to place an upper bound on Bg— Bg mix- 
ing. At this point we remark that our result corresponds 
to the Higgs sector specified above (supplemented by an 
additional 16 to avoid unwanted D-term breaking), which 
is minimal in its effect on the RG evolution of yt and leads 
to the weakest possible bound on Bg — Bg mixing. Chang- 
ing e.g. the last term in into a dimension-4 coupling 
involving a 126h Higgs field while keeping the low-energy 
parameters rrit and tan /3 fixed will strengthen our bound 
onBg — Bg mixing further. 

2 The top Yukawa coupling 

The large top Yukawa coupling yt suppresses the third- 
generation squark (and slepton) masses through its RG 
effects and generates the desired flavour-changing neu- 
tral bR-SR transitions. In both SO(IO) and SU(5) yt pos- 
sesses an infrared quasi-fixed point corresponding to a 
fixed point of the ratio yt/g, where g is the gauge coupling. 
From the values of mt and tan/?, which is the ratio of the 
VEVs of Hu and Hd, one can compute yt (x nit/ sin/3 at 
the electroweak scale and evolve it up to Mqvt, Mio and 
the fundamental scale near Mpi. The running above Miq 
is much stronger than between Mqut and Mio because 
the group-theoretical factors in SO(IO) are larger than in 
SU(5). Next we consider the critical Yukawa coupling , 
which corresponds to the quasi-fixed point in SO(IO), i.e. 
yt/g is constant above Miq. j/j is shown as a dashed curve 
in Fig. n with the two vertical lines indicating the scales 
il^GUT and MiQ. For yt < yt one consequently finds yt 
small at high energies. The solid line in Fig. ^ illustrates 
this situation for the following input parameters: 
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We have computed all relevant RG coefficients in SO(IO) 
and SU(5) using the general result of For the RGE's 
above Mqut we work in the leading logarithmic approx- 
imation, but include next-to-leading-order corrections in 
the MSSM. In this way we can account for electroweak 
threshold corrections, which become relevant if yt is close 
to yt - For instance, yt depends on the listed squark masses 
through these corrections. ai,2,s are the usual squared 
MSSM gauge couplings divided by An, 9t is the stop mix- 
ing angle and rrig is the gluino mass. One finds yt raised 
to yt for (mt,tan/3) — (180 GeV, 3) or for (m(,tan/3) = 
(184 GeV, 4) with the remaining parameters unchanged. 
For very small tan/3 < 2, the fixed-point values can be 
exceeded and the coupling typically becomes nonpertur- 
bative below the Planck scale. In this case the model loses 
its predictivity. Since further such low values of tan (3 are 
strongly disfavoured by LEP data, we require yt yt- 



3 Bg — Bs mixing 

The Bg — Bg mixing amplitude M12 can be expressed in 
terms of Wilson coefficients, which contain the short-dis- 
tance information, and matrix elements of local four-quark 
operators. The Standard Model contribution involves only 
a single operator: 

Ol = SLjf_,bL SL^bL- 

Its Wilson coefficient Gl is computed from the box di- 
agram involving two W bosons and top quarks. In the 
CMM model a new operator Or occurs, which is obtained 
from Ol by replacing the left-handed fields with right- 
handed ones. Since {Bg \Ol \ Bg) = (Bg \Or \ Bg), no new 
hadronic matrix elements are needed. Using the standard 
relativistic normalisation (S, IS,) — 2EV we can write 



A/12 



327r2 Mb, 



XUCL + CR){Bg\OL\Bg). (4) 



Here At = V'j* Vtb comprises the CKM elements. We define 
all Wilson coefficients and matrix elements at the reno- 
malisation scale fi — rrih, at which the Standard Model 
coefficient evaluates to Cl = 8.5 (5]. The leading contri- 
bution to Bg — Bg mixing in the CMM model arises from 
one-loop box diagrams with gluinos and squarks. The re- 
sult is 
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Here 



|^3| = |(7^3||f/r3| 



(5) 



(6) 



is the relevant combination of mixing matrix elements in 
the right-handed sdown sector, and 



[G{xi, xi) + G{xs,xs) - 2G{xi, xg)] 
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18 

[F{xj,,x~,) + Fixs,xs) -2F{x~,,x,)] (7) 



with Xq — rn1/ml and 5,6 denoting the right-handed 
squarks of the second and third generations. The functions 
F and G are defined in J5 . Note the twofold enhancement 
of Cr due to the large atmospheric mixing and the large 
strong coupling constant. This is, however, partially off- 
set by a smaller loop function 5^^^. The neutral Bg-meson 
mass difference is given by 



AMb, =2|Afi2| 



(8) 



The phase of M12 is responsible for mixing-induced CP 
violation. Note that the phase of A3 is undetermined, so 
that there is potentially large, but not predictable, CP 
violation in decay modes like Bg — > ipcf), tprj^''^ . A measure- 
ment of AMb, and the mixing-induced CP asymmetries 
in one of these decays will allow to determine both mag- 
nitude and phase of M12 and therefore of Cr. If the Bg 
oscillations are too rapid, these asymmetries cannot be 
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Fig. 1. RG evolution of yt- See text for explanation. 



measured. Then still the width difference in the sys- 
tem can be used to determine cos(argMi2) [7|. 

The input parameters for our analysis are j/t, rrig and 
the average squark mass rriq and the trilinear term aa of 
the first generation, all defined at the electroweak scale. 
They are evolved to the GUT scale and the SU(5) and 
SO(IO) GUT parameters are determined. All other param- 
eters are then determined from the RG flow back down 
to the electroweak scale. The flavour-changing effects in 
the CMM model grow with the sizes of the universal soft 
sfermion mass toq a-nd the universal trilinear term oq at 
the Planck scale. However, a larger mo also implies larger 
squark masses leading to a suppression of the gluino box 
function S^s'. The effect on Bs — Bs mixing is maximal 
for values of rriq around 800 GeV. Furthermore oo and 
mg, which for a given gluino mass determine the sfermion 
mass spectrum, are constrained by the lower bounds on 
these masses. We have further checked that no charge- 
or colour-breaking vacuum occurs. Finally the effect de- 
creases with increasing m^, so that the maximal effect is 
found for the experimental lower bound rrig — 195 GeV. 

The allowed area in the complex M12 plane for the 
CMM model is depicted in Fig. |2 For simplicity we ne- 
glect the 30% uncertainty from current lattice QCD de- 
terminations of {Bs \Ol\Bs)- Fixing the value of this ma- 
trix elements to the central value quoted in |B] results in a 
Standard Model prediction of 17.2 ps~^. The small black 
circle in Fig. |21 indicates this value and the large red disk 
denotes the range covered by the CMM model. The blue 
circle centered at the origin is the region excluded by the 
95 % CL experimental lower bound |S] 



AMb, > 14.4 ps" 



(9) 



We see that in total the mass difference can exceed its 
Standard-Model prediction by a factor of 16. 

In conclusion we have perfomed a RG analysis for the 
CMM model and computed the impact on Bg — Bg 
mixing. Our work complements and improves previous 
GUT-inspired analyses (see e.g. W), which supplement the 
MSSM with minimal flavour violation by bu-sn mixing at 
the electroweak scale. 
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Fig. 2. Allowed region in the complex 2Mi2 plane. The axis 
units are inverse picoseconds. 
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